arXiv:1504.08162vl [physics.optics] 30 Apr 2015 


Synthetic gauge fields for light beams in optical resonators 


Stefano Longhi 

Dipartimento di Fisica, Politecnico di Milano and Istituto di Fotonica e Nanotecnologie del Consiglio Nazionale delle 
Ricerche, Piazza L. da Vinci 32, F20133 Milano, Italy (stefano.longhi@polimi.it) 


Compiled May 1, 2015 

A method to realize artificial magnetic fields for light waves trapped in passive optical cavities with anamorphic 
optical elements is theoretically proposed. In particular, when a homogeneous magnetic field is realized, a 
highly-degenerate Landau level structure for the frequency spectrum of the transverse resonator modes is 
obtained, corresponding to a cyclotron motion of the optical cavity field. This can be probed by transient 
excitation of the passive optical resonator. © 2015 Optical Society of America 
OCIS codes: 140.4780, 070.5753, 000.1600 


Introduction. Realizing magnetic effects with photons at 
small spatial scales is a challenge. However, it would 
be extremely useful from both fundamental and ap¬ 
plied aspects, enabling for instance the possibility to 
study quantum Hall effects for light, to built compact 
integrated optical isolators and to realize robust light 
transport via topologically-protected states [1-6]. One 
possibility relies on the of synthesis of artificial gauge 
fields, which mimic magnetic or pseudo-magnetic fields 
for photons. This has been suggested and experimentally 
demonstrated in a series of recent works (see [6-19] and 
references therein) using twisted [7,8] or strained [9,10] 
waveguide lattices, coupled-resonator optical waveguide 
structures [6,11-13], photonic lattices with suitable tem¬ 
poral modulation of the lattice parameters [14-17], and 
driven dissipative lattices of polaritons [18,19]. Such pre¬ 
vious studies have been mainly concerned with tight- 
binding lattice structures in either square or honeycomb 
geometries, where artificial gauge fields are generally 
synthesized by introduction of effective Peierls phases in 
the hopping rates among adjacent lattice sites. Another 
class of optical systems where light waves can behave 
like quantum particles in external potentials is provided 
by optical resonators [20], where the transverse mode 
spectrum and corresponding beam dynamics can be con¬ 
trolled by proper design of the optical potential [21-25]. 

In this Letter we suggest a method to realize artificial 
magnetic fields for light waves trapped in passive optical 
cavities with anamorphic optical elements. In particular, 
for a homogeneous magnetic force a highly-degenerate 
Landau level structure for the frequency spectrum of the 
transverse resonator modes is obtained. 

Quantum particle in a magnetic field and Landau lev¬ 
els. The two-dimensional motion of a quantum parti¬ 
cle of mass m and charge q, constrained on the plane 
(x, y) in the presence of external magnetic and electric 
fields B = V X A, E = —VH described by the vectorial 
and scalar potentials A = Ax{x,y)ux + Ay{x,y)uy and 
V{x,y), is governed by the Schrodinger equation 

= T [{p^ - qA^f + {Py - qAy)^] ^+qV'ip = 

( 1 ) 


for the particle wave function fj{x,y,t), where px = 
—ihdx and py = —ifidy are the canonical momenta. In 
particular, for a uniform magnetic field (B = 

E = 0) in the Landau gauge one can assume Ax = 
—Boy, Ay = 0, V = 0, leading to the Hamiltonian 
PL = (l/2m) [{px — qAxY well-known, the 

eigenfunctions of such an Hamiltonian are a product of 
momentum eigenstates in the x direction and harmonic 
oscillator eigenstates in the y direction, namely 


'Ip = 'ipn= exp{ikxx)Hn ~ 2/o) j ^ 

X exp [-mu)c{y - yof /2h] , (2) 

where ujc = qB^jm is the cyclotron frequency of clas¬ 
sical motion and yo = fikxj [rnu^f). The corresponding 
energies En, 'H'lpn = form a set of Landau levels, 

En = hujc{n -h 1/2), which do not depend on the quan¬ 
tum momentum kx and are thus highly degenerate levels. 
Here we show that an optical Schrodinger equation like 
Eq.(l) arises in astigmatic optical cavities. 

Synthetic magnetic field in an optical resonator. Let 
us consider the Fabry-Perot optical resonator shown in 
Fig. 1(a). It comprises two high-reflectivity mirrors, a 
plane mirror (at the left side) and an aspherical (quasi- 
prismatic) mirror at the right side. The aspherical mirror 
has a reflection coefficient r{x,y) = ex.p[ikxG{y)], where 
the function G{y) is left arbitrary at this stage. Note 
that this is a quasi-prismatic mirror, i.e. it has a strictly 
prismatic shape in the x direction and an arbitrary struc¬ 
ture in the y direction. Four lenses are placed between 
the two mirrors, namely two cylindrical lenses (LI) with 
focal length / and two spherical-cylindrical lenses (L2) 
with focal lengths fx = f and fy = f /2 m the x and y di¬ 
rections, respectively. The distances between the various 
adjacent lenses and mirrors are indicated in the figure. 
A thin plate P, acting as a phase filter, is also placed 
close to one of the two lenses LI, as depicted in Fig. 1(a), 
with a transmission function t{x,y) = exp[i(^(x, ^)]. A 
monochromatic optical field E{x,y,z,t) at frequency 
uu = kc = 2i\cj\ trapped in the resonator is described 
by the superposition of counter-propagating waves as 
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Fig. 1. (Color online) (a) Schematic of the optical resonator 
with astigmatic optical elements that realizes a synthetic 
magnetic held in the transverse plane (LI: cylindrical lens; 
L2: spherical-cylindrical lens), (b) Round-trip propagation of 
the intracavity counter-propagating envelopes F and B. 

E = F(x, z) exp(ici;t — ikz) + B{x^ y, z) exp(ici;t + ikz). 
The propagation of the envelopes F and B at different 
planes inside the cavity can be readily obtained by ap¬ 
plication of the generalized Huygens-Fresnel integral for 
astigmatic beams [20,26]. In particular, each of the two 
spherical-cylindrical lenses L2 realize a one-dimensional 
Fourier transformation in the x direction and an imaging 
step in the y direction between the two planes far apart 
by / from each lens. Using such a property, it can be 
readily shown that the following relations hold between 
the forward and backward propagating envelopes Fi ^2 
and Bi ^2 at the two planes 1 and 2 shown in Fig. 1(b) 

B2{x,y) =V:cF2{x + fG{y),y) , Fi{x,y) = fo,Bi{x,y). 

(3) 

where the symmetry operator Vx is dehned by 
Vx'^ix^y) = '0(—x,^). In writing Eq.(5), we assumed 
100 % mirror reflectance and no injected signal into the 
resonator; the effects of output coupling and signal injec¬ 
tion will be discussed below. Taking into account diffrac¬ 
tion and phase effects between the two planes 1 and 2, 
one can also write 

F 2 (x, y) = t{x, y)VFi (x, y) , Bi (x, y) = Vt{x, y)B 2 (x, y) 

( 4 ) 

where V = exp[—'i//(2/c)Vf] is the diffraction operator 
for a propagation length /, and Vf = 9^ + dy is the 
transverse Laplacian. From Eqs.(3) and (4) the resonator 
round-trip propagator /C at the plane 2, which connects 
the amplitude F 2 {x^y) of forward waves at successive 
transits in the cavity can be readily derived and reads 

t = tVVxBtVxS (5) 

where we have further introduced a shift operator 5, 
defined by S 2 p{x,y) = 2 p{x + fG{y),y). Provided that a 
single longitudinal frequency of the cavity is excited and 


assuming the reference wave number k to coincide with a 
longitudinal cavity wave number, the beam dynamics at 
successive transits inside the resonator can be described 
by the following map 

{x, y) = /£</)(") {x, y) (6) 

where (a;, y) = FF{x,y) is the envelope of the pro- 
gressive wave at plane 2 and at the n-th transit in the 
cavity. In the limit / ^ 0, (/9 ^ 0, and G 0, i.e 
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Eig. 2. (Color online) Cyclotron motion of a Gaussian beam 
in the resonator under a homogeneous magnetic field: (a) 
Snapshots of the transverse intensity distribution at plane 2 
at a few successive round trip numbers n. (b) and (c): tra¬ 
jectory of the beam center of mass (x, y) in the transverse 
plane. Parameter values are given in the text. 



Eig. 3. (Color online) Cyclotron motion of a Gaussian beam 
under a gradient magnetic held, (a) and (b) show the trajec¬ 
tory of the beam center of mass (x, y) in the transverse plane. 
Parameter values are given in the text. 

provided that the diffractive and the phase effects (in¬ 
troduced by the free-space propagation for a length I 
and by the phase plate and prismatic mirror) are weak, 
the resonator round-trip operator 1C reduces at leading 
order to the identity operator X, because VxBx = X. In 
such a case, the held envelope (x, y) undergoes a slow 
change after each round trip and one can set fC 
where the leading order correction iCl to the identity 
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operator is obtained from Eq.(5) by first-order power 
expansion of the operators P, t and S. After continua¬ 
tion of the round-trip index n and setting (x, y) — 

(j)^'^\x^y) {d(j)/dn) with (j)[x^y^n) = (j)^'^\x^y)^ one 

then obtains the evolution equation (mean-field equa¬ 
tion) i{d(j)/dn) = which can be cast in the follow¬ 

ing form 


ih-p— = -— I —iTi— - qA 


dn 2m 


dx 




In Eq.(7) we have set ^{x^y^n) = (j)*{x,y,n) and 


n 

m 

qAx 

qV{x,y) 


1/k = A/(27r) 

1 /( 2 /) 

fG{y) 

21 

--[(p{x,y) +ip{-x,y)] 


( 8 ) 

(9) 


( 10 ) 


PG^y) 

4/ 


( 11 ) 


The mean-field equation (7) clearly corresponds to the 
Schrodinger equation (1) for a quantum particle in an 
external electric and magnetic field with = 0, pro¬ 
vided that the formal substitutions (8-11) are made in 
the Schrodinger equation. Note that the round-trip num¬ 
ber n in Eq.(7) corresponds to the time variable normal¬ 
ized to the cavity round-trip time Tr = 2I/e/c, where Lg 
is the optical resonator length. Note also that the mag¬ 
netic field is controlled by the phase profile G{y) of the 
quasi-prismatic mirror, and turns out to be generally in¬ 
homogeneous in the y direction. The scalar potential V 
is in turn controlled by both the phase profile G{y) and 
the phase profile ip{x,y) of the phase plate P according 
to Eq.(ll). In particular, with the choice 


(p{x,y) = -kfG‘^{y)/{8l) (12) 


the potential V vanishes, and Eq.(7) describes the mo¬ 
tion of a quantum particle in an inhomogeneous mag¬ 
netic field. Note that ip defined by Eq.(12) does not de¬ 
pend on X and it thus describes the effect of a (generally 
aspherical) cylindrical lens. 

Transverse mode speetrum and beam dynamies. Let us 
first consider the case of a homogeneous magnetic field. 
This is obtained by assuming a linear gradient for G, 
i.e. G{y) = ay^ which corresponds to a magnetic field 
qB^ = af /{21). Note that in this case a quadratic phase 
if [Eq.(12)] is required to cancel the scalar potential E, 
i.e. the phase plate P in Eig.l(a) can be simply replaced 
by a converging cylindrical lens with focal length /' in 
the y direction given by /' = 2l/{fa)‘^. The transverse 
mode spectrum corresponds to a set of highly-degenerate 
Landau levels at frequencies = cj -h 27rAz/(l/2 -h n) 
with frequency separation 


Az/ = 


UJc 

2'kTr 


27tTr 


af 


( 13 ) 


where n = 0,1,2,... is the Landau level index and 
Auax = 1/Tr = c/(2I/e) is the frequency separation 


of the cavity axial modes. The corresponding transverse 
modes are given by Eq.(2). Note that the Landau degen¬ 
eracy of transverse modes is rather distinct than degen¬ 
eracy found in the spectrum of other broad area opti¬ 
cal cavities (such as plane-plane or confocal resonators) 
generally employed to study transverse patterns in lasers 
and nonlinear optical systems [27,28]. A distinctive fea¬ 
ture here is the possibility to observe a cyclotron mo¬ 
tion as a result of mode degeneracy. This is shown, as 
an example, in the numerical simulations of Eig.2 for 
resonator parameters / = 4 cm, / = 0.4 cm, corre¬ 
sponding to a total cavity length Lg — 4/ -h / — 16.4 
cm, probed at the wavelength A = 633 nm. Assuming 
G{y) = ay with a = 1 m“^, the frequency separa¬ 
tion of the set of degenerate Landau modes is given by 
Az/ 5.82 MHz, corresponding to a cyclotron motion 
that occurs after rig = 27r/ {af) 157 round trips in the 
resonator, i.e. with a period Tg = UcTr = l/Az/ 172 
ns. Eigure 2(a) shows snapshots of the transverse beam 
intensity evolution at successive round trips as obtained 
by numerical solution of the map (6), taking as an ini¬ 
tial condition an off-axis Gaussian field distribution with 
spot size w and transverse displacements xi and ^i, 
namely (x exp[—(x — x\p jvf^ — {y — yi)‘^/w^] with 
w = 300 /im, xi = 0.5 mm and = 1.5 mm. The trajec¬ 
tory of the beam center of mass in the transverse (x, y) 
plane, taken every 5 round trips, is depicted in Eig.2(b) 
and clearly corresponds to a nearly circular path of ra¬ 
dius R yi, according to a semiclassical analysis of the 
Hamiltonian R with the chosen Landau gauge. Interest¬ 
ingly, for a weak inhomogeneous magnetic field in the 
y direction, corresponding to {d?G/dy‘^) ^ 0, the beam 
path is given by the superposition of a relatively fast cir¬ 
cular motion around a point (the guiding center) and a 
relatively slow drift of this point along the x-axis. The 
weak drift of the guiding center is shown in Eig.3 for a 
quadratic phase profile G{y) = ay — ^y^ with a = 1 m“^ 
and P = 30 m“^; other parameters are as in Eig.2. 

The cyclotron motion depicted in Eigs.2 and 3 could be 
experimentally observed by monitoring the intensity dis¬ 
tribution in a resonator driven by an external beam. Indi¬ 
cating by T ^ I the transmittance of the coupling plane 
mirror [at the left side in Eig.l(a)] and by E\^{x^y) the 
injected field amplitude at a frequency in resonance with 
one axial mode of the resonator and slowly varying over 
one round trip, the evolution of the intracavity field is 
now described by the following modified map [compare 
with Eq.(6)] 


^(”+1) (t) = + VteIP - . (14) 

Similarly, the mean-field equation (7) takes the form 



T 

ih—tp + iUVrEin. 


2m 


+ Vip 


(15) 


As an example, Eig.4 shows the evolution of the beam po- 
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sition inside the resonator when a Gaussian beam is cou¬ 
pled into the cavity via the plane mirror (T = 1% trans¬ 
mittance) and slowly switched off. In the simulations, we 
assumed = f{n) exp[—(x —— {y — yiY/w‘^]^ 
where the function f{n) = (1/2)[1 — tanh((n — no)/An)] 
describes the switch off of the the injected signal at the 
round trip no [see the dashed curve in Fig.4(b)]. Pa¬ 
rameter values used in the simulations are the same as 
those in Fig.2, except for a = 2 m“^. For n < no, the 
intracavity field is stationary and it shows a character¬ 
istic annular shape [Fig.4(a)], which is the signature of 
a cyclotron dynamics: the annular shape at steady-state 
operation arises from a balance between circular drift 
induced by the synthetic gauge field and continuous sig¬ 
nal feeding at x = xi, ^ = ^i. As the external field is 
switched off (n > no) and signal feeding vanishes, the 
intracavity field decays in time undergoing a cyclotron 
motion, see Fig.4(b). Note that the cyclotron period is 
l/An 86 ns, which is smaller than the cavity decay 
rate Tr/T 109.4 ns. Hence the spiraling of the decay¬ 
ing output light beam when the holding beam is switched 
off should be observable by time-resolved transverse pat¬ 
tern evolution measurements [29,30]. 



y (mm) round trip number n 


Fig. 4. (Color online) Cyclotron motion in a resonator with 
injected Gaussian beam, (a) Stationary intensity distribution 
at plane 2 under signal injection (n < no), (b) Trajectory 
of the beam center of mass (x, y) when the injected field is 
switched off. The behavior of the intracavity power (in arbi¬ 
trary units) is also shown, together with the decaying ampli¬ 
tude/(n) = (1/2)[1—tanh((n —no)/An)] of the injected field 
(dashed curve) with no = 200 and An =10. Other parameter 
values are given in the text. 

Conclusions. Artificial magnetic fields for photons can 
be realized in the transverse plane of an optical res¬ 
onators with astigmatic elements. Such a result sheds 
new light into the rapidly growing field of synthetic mag¬ 
netic fields in optical systems and might suggest further 
investigations. For example, the degenerate Landau level 
structure of the resonator can provide a platform for the 
study of laser mode competition and transverse pattern 
formation [27,28] in presence of a circular drift (cyclotron 
motion). Also, by including a periodic phase mask in 
the cavity, the coherently-driven astigmatic resonator 
could implement a driven-dissipative two-dimensional 
photonic lattice with a synthetic gauge field [19], thus 
providing a viable route toward the experimental obser¬ 
vation of optical analogues of the anomalous and the 
(integer) quantum Hall effect. 
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